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We study reactions between neutron-rich 132 Sn nucleus and 96 Zr within a dynamic microscopic 
theory at energies in the vicinity of the ion-ion potential barrier peak, and we compare the properties 
to those of the stable system 124 Sn+ 96 Zr. The calculations are carried out on a three-dimensional 
lattice using the density-constrained Time-Dependent Hartree-Fock method. In particular, we cal- 
culate the dynamic excitation energy E*(t) and the quadrupole moment of the dinuclear system, 
Q2o(t), during the initial stages of the heavy-ion collision. Capture cross sections for the two re- 
actions are analyzed in terms of dynamic effects and a comparison with recently measured data is 
given. 

PACS numbers: 21.60.-n,21.60.Jz 
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I. INTRODUCTION 



Heavy-ion reactions at radioactive ion beam (RIB) fa- 
cilities allow us to form new exotic neutron-rich nuclei 
and to study their physical properties. Examples include 
experiments with neutron-rich 132 Sn beams on targets of 
64 Ni [l], @] and of 96 Zr Q. Another experimental frontier 
is the synthesis of superheavy nuclei in cold fusion re- 
actions involving spherical closed-shell 208 Pb targets [4j 
and in hot fusion reactions with deformed actinide nu- 
clei J5J. These experiments present numerous challenges 
for a theoretical description, in particular for dynamic 
microscopic theories. 

At relatively large impact parameters, heavy-ion re- 
actions are dominated by deep inelastic collisions in 
which the nuclei make only brief contact. The reac- 
tion products have mass and charge similar to projec- 
tile and target, but the energy may be strongly damped. 
At smaller impact parameters, an intermediate dinuclear 
system is formed. Entrance-channel heavy-ion poten- 
tials have been calculated in various models, including 
the macroscopic-microscopic method with five shape pa- 
rameters [6|, and the energy density functional method 
with extended Thomas- Fermi approximation [7]. Fur- 
thermore, dynamical models [8HlC( show that if the din- 
uclear system is able to move inside the saddle point, 
capture occurs which may lead to the formation of a 
compound nucleus with compact shape. During capture 
the energy of relative motion of the ions gets converted 
into intrinsic excitation energy E*(t). If a compound 
nucleus is formed it will subsequently decay either by 
particle evaporation or by fission. On the other hand, if 
the dinuclear system separates before crossing the saddle 
point, the reaction process is called quasifission. In this 
case, the shape of the dinuclear system is very elongated 
(large quadrupole moment). Experimentally, it is pos- 
sible to separate fusion-fission from quasifission by mea- 
suring the angular anisotropy of the fragments. In the 



collision of very heavy ions such as 132 Sn+ 96 Zr consid- 
ered in this paper, studies of fusion reactions are com- 
plicated by the competition with quasifission and fusion- 
fission events which hinder the formation of evaporation 
residues. 

The time-dependent Hartree-Fock (TDHF) theory pro- 
vides a useful foundation for a fully microscopic many- 
body theory of heavy- ion collisions in the vicinity of the 
Coulomb barrier [ll|, [l2| ■ Partly because of the recent 
breakthroughs in microprocessor technology, it has be- 
come feasible to perform TDHF calculations on a three- 
dimensional (3D) Cartesian grid with no symmetry re- 
strictions and with much more accurate numerical meth- 
ods. At the same time the quality of effective interactions 
has also been substantially improved [13, [Tj] . These de- 
velopments allow for the testing of the time-dependent 
mean-field approach to nuclear reactions without any nu- 
merical uncertainties |15l - tla |. The TDHF code used in 
these calculations utilizes the full Skyrmc interaction, in- 
cluding all of the time-odd terms in the mean field Hamil- 
tonian [UGl]. 

During the past several years, we have developed the 
density-constrained (DC) TDHF method (DC-TDHF) 
for calculating heavy- ion potentials [20J. We have ap- 
plied this method to calculate fusion cross sections above 
and below the barrier for a number of systems: The first 
application was for the 132 Sn+ 64 Ni [HHl] system. The 
fusion cross section at the lowest projectile energy has 
been re-measured [2[ and now agrees remarkably well 
with our calculations. We have also performed calcu- 
lations for 64 Ni+ 64 Ni |3 and for 16 O+ 208 Pb[24]. In all 
these cases, we have found very good agreement between 
the measured fusion cross sections and the DC-TDHF 
results. Very recently, we have carried out a microscopic 
dynamical study of the astrophysical triple-a reaction to 
form a resonant state of 12 C [25[ and a study similar to 
the one presented here for reactions involving superheavy 
formations [26| , using the same approach. 

In the present paper, we study reactions between 



the neutron-rich 132 Sn nucleus and 96 Zr at energies in 
the vicinity of the ion-ion potential barrier peak, and 
we compare observables to those of the stable system 
124 Sn+ 96 Zr. The dynamic microscopic calculations are 
carried out on a 3D Cartesian lattice using both unre- 
stricted TDHF and DC-TDHF methods. This is by far 
the heaviest neutron-rich system we have investigated so 
far, and the microscopic numerical calculations with the 
added DC are computationally very intensive. 

This paper is organized as follows: in Section II we 
summarize the Formalism (DC-TDHF, dynamic excita- 
tion energy, capture cross section) . In Section III numer- 
ical results are presented. In particular, we show contour 
plots of the mass density of the dinuclear system and dis- 
cuss the dynamic quadrupole moment Q2o(t) during the 
initial stages of the collision. We also calculate the heavy- 
ion interaction potential V(R) and demonstrate that in 
these very heavy systems the barrier height and width 
increase dramatically with increasing beam energy. In- 
teraction barrier heights and positions are also deduced 
from unrestricted TDHF runs. We examine the dynamic 
excitation energy E*(t) during the initial stages of the 
collision and compare it to the excitation energy of the 
compound nucleus in its ground state, E* = E c , m ,+Q gg . 
Finally, capture cross sections for the two reactions are 
analyzed in terms of dynamic effects, and a comparison 
with recently measured capture-fission data [3J is given. 
The conclusions are presented in Section IV. 



II. FORMALISM: DC-TDHF METHOD, 
DYNAMIC EXCITATION ENERGY, CAPTURE 

Recently, we have developed a method to extract ion- 
ion interaction potentials directly from the TDHF time- 
evolution of the nuclear system. In our DC-TDHF ap- 
proach [20j . the TDHF time-evolution takes place with 
no restrictions. At certain times during the evolution the 
instantaneous density is used to perform a static Hartree- 
Fock minimization while holding the total proton and 
neutron density of the dinuclear system constrained to 
be the instantaneous TDHF density. This provides us 
with the TDHF dynamical path in relation to the multi- 
dimensional static energy surface of the combined nuclear 
system. In the DC-TDHF method the ion-ion interaction 
potential is given by 



V(R)=E DC (R)-E Al -E A , 2 



(1) 



where -Edc is the DC energy at the instantaneous separa- 
tion R(t), while -Eai and E\ 2 are the binding energies of 
the two nuclei obtained with the same effective interac- 
tion. The interaction potentials calculated with the DC- 
TDHF method incorporate all of the dynamical entrance 
channel effects such as neck formation, particle exchange, 
internal excitations, and deformation effects. While the 
outer part of the potential barrier is largely determined 
by the entrance channel properties, the inner part of the 



potential barrier is strongly sensitive to dynamical effects 
such as particle transfer and neck formation. 

For the calculation of the ion- ion separation distance R 
we use a hybrid method as described in Ref. [24| . At large 
distances where a still visible neck allows us to identify 
two fragments we compute it as distance of the center of 
mass of the ions. For more compact configurations, we 
compute R from the mass quadrupole moment Q20 as 
R = r o\/|Q20 where r is a scale factor to connect the 
definition smoothly to the large-distance region. 

In heavy-ion reactions, the total capture cross section 
consists of the following terms [3] : 



-'capt 



o'er + ctqf + off 



(2) 



where oer, cqf, off denote the evaporation residue cross 
section, the quasifission cross section, and the fusion- 
fission cross section. In the reaction of light and medium- 
heavy ions, the fission barriers of the pre-compound sys- 
tem are so high that fission contributions are negligible. 
In this case we have <7 cap t ~ o'er- On the other hand, 
for the reaction of massive nuclei like 132 Sn+ 96 Zr, the 
pre-compound system is an excited state of the actinidc 
nucleus 228 Th with a fission barrier of only about 6 MeV; 
thus we expect sizable fission contributions, and the evap- 
oration residue cross section is expected be be rather 
small. The number of quasifission events increases dra- 
matically with the product of the charge numbers Z\Z% 
and with the orbital angular momentum I in the entrance 
channel. Another reason for the decreasing yield of ER 
formation is that a heated and rotating CN may fission 
(fusion-fission) . 

Ion-ion interaction potentials calculated using DC- 
TDHF correspond to the configuration attained during 
a particular TDHF collision. As mentioned above, for 
light and medium mass systems as well as heavier sys- 
tems for which fusion is the dominant reaction channel, 
DC-TDHF gives the fusion barrier with an appreciable 
but relatively small energy dependence. On the other 
hand, for reactions involving massive systems fusion is 
not the dominant channel at barrier top energies. In- 
stead the system sticks in some dinuclear configuration 
with possible break-up after exchanging a few nucleons. 
The long-time evolution to break-up is beyond the scope 
of TDHF due to the absence of quantum decay processes 
and transitions. As we increase the energy above the bar- 
rier this phenomenon gradually changes to the formation 
of a truly composite object. This is somewhat similar 
to the extrapush phenomenon discussed in phenomeno- 
logical models. For this reason the energy dependence 
of the DC-TDHF ion-ion potential barriers for these sys- 
tems is not just due to the dynamical effects for the same 
final configuration but actually represents different final 
configurations. 

Theoretically, the calculation of the total capture cross 
section is similar to the calculation of the fusion cross 
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with the understanding that the ion-ion interaction po- 
tential used in the calculations distinguishes the two 
events. In practice, the potential barrier penetra- 
bilities Tl at E c . m . energies below and above the 
barrier are obtained by numerical integration of the 
Schrodinger equation for the relative coordinate R using 
the well-established Incoming Wave Boundary Condition 
(IWBC) method (HHi]. This Schrodinger equation con- 
tains the heavy-ion potential V(R) given in Eq. (TTJ and 
the centrifugal potential. In the IWBC calculations the 
summation over L in Eq. [3] is continued until the contri- 
bution becomes negligible to the total cross section. As 
we shall discuss below, we can also determine the max- 
imum value of L by performing TDHF calculations for 
non-zero impact parameters. The optimal way to study 
the problem would be to perform DC-TDHF calculations 
for different L values, however for heavy systems the com- 
putational cost for doing this is very large. Finally, for 
the calculation of capture cross sections it is possible to 
use a coordinate dependent effective mass p(R) as de- 
scribed in Ref. |24|. The effect of using a coordinate 
dependent mass is to modify the inner part of the ion- 
ion potential, particularly at low subbarrier energies. For 
the energies studied here we have found this effect to be 
very small for capture cross sections. 

Taking up the strategy proposed in [30(, we have also 
developed a new microscopic approach for calculating dy- 
namic excitation energies E* (t) of systems formed during 
heavy-ion collisions [31|. For this purpose, we divide the 
conserved TDHF energy into a collective and intrinsic 
part, and we assume that the collective part is primarily 
determined by the density p(r, t) and the current j(r, t). 
Consequently, the excitation energy can be written in the 
form 
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FIG. 1. (color online) TDHF calculations for 132 Sn+ 96 Zr. 
Contour plots of the mass density at the distance of closest 
approach in a central collision, calculated at three different 
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energies. 



E*(t) = £ T dhf - E con (p(t),j(t)) 



(4) 



where £tdhf is the total energy of the dynamical system, 
which is a conserved quantity, and -Ecoii represents the 
collective energy of the system. The collective energy 
consists of two parts 

-Econ (t) = Skin (p(t), j(t)) + E DC (p(t)) , (5) 

where -Ekin represents the kinetic part and is given by 



Skin (p(t), j(t)) = 



d 3 rf(t)/p(t), 



(6) 



which is asymptotically equivalent to the kinetic energy 
of the relative motion, ^pR 2 , where p is the reduced 
mass and R(t) is the ion-ion separation distance. The en- 
ergy -Edc is the lowest-energy state of all possible TDHF 
states with the same density and is required to have zero 



excitation energy. The dynamics of the ion-ion separa- 
tion R(t) can be extracted from an unrestricted TDHF 
run. Using E*(t) and R(t), we can deduce the excitation 
energy as a function of the distance parameter, E*(R). 



III. RESULTS 

The numerical calculations are carried out on a 
3D Cartesian lattice using the Basis-Spline collocation 
method to represent derivative operators with high accu- 
racy. For the 132 ' 124 Sn+ 96 Zr reactions studied here, the 
lattice spans 50 fm along the collision axis and 30 — 42 fm 
in the other two directions, depending on the impact pa- 
rameter. The lattice spacing is 1.0 fm in all directions. 
We utilize the full Skyrme interaction (SLy4) [l3[ includ- 
ing all of the time-odd terms in the mean field Hamilto- 
nian llal. without the cm. correction as described in 



Ref. [32| . The two nuclei are placed at an initial sepa- 
ration of 22 fm. First we generate highly accurate static 
HF wave functions for the two nuclei on the lattice, which 
are then boosted and time-propagated with a time step 
At = 0.4 fm/c. The computation of the dynamic exci- 
tation energy and the heavy-ion potential is numerically 
very intensive, primarily due to the DC calculation. In 
a typical DC-TDHF run, we utilize a few thousand time 
steps, and the DC is applied every 20 time steps. To dis- 
tinguish between deep-inelastic and capture reactions, we 
have also performed several unrestricted TDHF runs for 
the 132 Sn+ 96 Zr system above the barrier. The numerical 
accuracy of the static binding energies and the deviation 
from the point Coulomb energy in the initial state of the 
collision dynamics is on the order of 50 — 200 keV. The 
accuracy of the DC calculations is commensurate with 
the accuracy of the static calculations. 
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A. Dynamic quadrupole moment, Interaction 
Barrier, and Capture Barrier 

In Fig. Q] we show contour plots of the mass density 
at the distance of closest approach in a central collision. 
These density plots reveal that at E c . m . = 195 MeV the 
nuclear surfaces barely touch; this energy corresponds to 
the interaction barrier. At E c . m . = 210 MeV we still 
notice a density configuration with two separate cores. 
Only at energies E c . m . = 230 MeV and above, a single- 
core composite system emerges, albeit with very large 
elongation. The large elongation of the composite sys- 
tem is readily apparent if one plots the intrinsic mass 
quadrupole moment 



O20 (*) = \f^-Jd\p(v, t)(2z 2 - x 2 - y 2 



(7) 



as a function of time, see Fig.[2j Also shown, for compar- 
ison, is the static quadrupole moment of the compound 
nucleus 228 Th which is more than three times smaller. 
Furthermore, the plot shows that central collisions at en- 
ergies -Ec.m. > 230 MeV lead to capture while the nuclei 
bounce back from each other at Ec.m. = 220 MeV and 
below (deep-inelastic collision). 

As we have discussed earlier, for systems leading to su- 
perheavy formation the evaporation residue cross section 
is customarily represented in terms of the various phases 
of the reaction process as 



O'ER — "capture 
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survival ; 



(8) 



where o'er denotes the evaporation-residue cross section 
for the superheavy system, cr C apture is the capture cross 
section for the two-ions, Pcn is the probability of form- 
ing a compound nucleus, and Survival is the probability 
that this compound system survives various breakup and 
fission events. The calculations presented here can only 
address the capture cross section for these systems since 
the subsequent reaction possibilities are beyond the scope 



t(10" 21 s) 



FIG. 2. (color online) Intrinsic mass quadrupole moment as 
a function of time. 



of the TDHF theory. For most light systems, for which 
fusion is the dominant reaction result, er C apture and oer 
are essentially the same and equal to the fusion cross 
section, Of US i on . For reactions involving superheavy for- 
mations we instead have 



''capture — O-QF + ^FF + O'ER , 



(9) 



where ctqf and off denote the quasi-fission and fusion- 
fission cross sections, respectively. For these reactions 
the evaporation residue cross section, o'er, is very small 
and therefore the capture cross section is to a large ex- 
tent equal to the sum of the two fission cross sections. 
Furthermore, the distinction between deep-inelastic re- 
actions and quasi-fission is somewhat difficult and usu- 
ally achieved by setting windows for fragment masses of 
Af = ^4cn/2 ± 20 and on their kinetic energy. 

In Fig.[3]we show the ion-ion potential V(R) for a cen- 
tral collision of 132 Sn+ 96 Zr, calculated at four different 
E c .m. energies using the DC-TDHF method. The dotted 
part of the potential line calculated at E c . m . = 220 MeV 
corresponds to the outgoing trajectory (the nuclei bounce 
back in a deep-inelastic collision). Our results demon- 
strate that in these very heavy systems the barrier height 
and width increase dramatically with increasing energy 
E c .m.- In fact, at higher energies the potential becomes 
almost flat. This is the first n-rich system in which we 
have observed this behavior. By contrast, DC-TDHF 
calculations for light ion systems such as 16 0+ 16 show 
almost no energy-dependence even if we increase E cm . 
by a factor of four 3l| . Even in reactions between a light 
and a very heavy nucleus such as 16 O+ 208 Pb, we see only 
a relatively small energy dependence of the barrier height 
and width [24| . For comparison, we have also plotted the 
phenomenological double- folding potential [33|, |34[ which 



uses the ground state densities of the two nuclei and keeps 
them frozen. This potential is energy-independent and 
has been calculated using the M3Y effective NN inter- 
action [39] and static 2-D HFB densities. We observe 
that the double-folding potential yields a potential bar- 
rier which is fairly similar to that of the DC-TDHF poten- 
tial at E c . m . = 230 MeV; however, because of the frozen 
density approximation, the potential exhibits too much 
attraction at smaller distances. Another difference is that 
the ground state densities in the double folding method 
correspond to well bound nucleons. The TDHF densities, 
on the other hand, have much wider tails because they 
cover dynamically excited nucleons. This puts some part 
of the excitation energy into the DC-TDHF energy and 
hence into the ion-ion potential. 
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FIG. 3. (color online) DC-TDHF calculations for the neutron- 
rich system 132 Sn+ 96 Zr. The potential barriers V(R) at four 
Ec.m. energies are obtained using Eq. |T}. Also shown is 
the point Coulomb potential. The phenomenological double- 
folding potential (dashed black curve) is given for comparison. 

In Fig. [J] we show the ion-ion potential barriers in de- 
tail. A comparison is made between the neutron-rich sys- 
tem 132 Sn+ 96 Zr and the stable system 124 Sn+ 96 Zr. We 
find that the potential barriers of the neutron-rich sys- 
tem are systematically 1 — 2 MeV higher than those of 
the stable system. We emphasize again that only the po- 
tential barriers calculated at energies E c . m . > 230 MeV 
lead to a true composite system with overlapping cores 
(capture) while the potential barriers calculated at ener- 
gies -E c .m. < 230 MeV correspond to a dinuclear system 
where both nuclei maintain separate cores (deep- inelastic 
collisions) . 

The interaction barrier Vj is defined as the energy to 
bring the two colliding nuclei into contact J3J. This en- 
ergy can be inferred from contour plots of the TDHF 
mass density at the distance of closest approach in a cen- 
tral collision, see Fig.[TJ For 132 Sn+ 96 Zr, we find an inter- 
action barrier height Vj = Ec.m. = 195 MeV, and the cor- 



responding distance of closest approach Rj = 14.77 fm. 
In Table U we summarize the interaction barriers and ion- 
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FIG. 4. (color online) Comparison of the heavy-ion barriers 
for the neutron-rich system 132 Sn+ 96 Zr (solid lines) and the 
stable system 124 Sn+ 96 Zr (dashed lines). The potential bar- 
riers are obtained with the DC-TDHF method at five _E c . m . 
energies. 



ion potential barriers (capture barriers) for the two sys- 
tems and their corresponding positions in i?-space. While 
the DC-TDHF barriers are fairly similar, we observe large 
differences (9 MeV) in the interaction barriers of the two 
systems: the additional neutrons in 132 Sn give rise to 
a larger attractive potential which causes the nuclei to 
snap together at lower energy. For TDHF collisions of 



TABLE I. Interaction barrier heights Vi (energy to bring the 
two colliding nuclei into contact) and barrier positions Ri 
calculated with unrestricted TDHF at zero impact parameter. 
Also given are the energy-dependent ion-ion potential barrier 
heights Vb and positions Rb determined with the DC-TDHF 
method. 



Reaction 


Vi (MeV) 


Ri (fm) 


V B (MeV) 


Rb (fm) 


132 Sn+ 96 Zr 


195 


14.77 


211.4^ 

215. 0~ 


13.03^. 
12.56~ 


124 Sn+ 96 Zr 


204 


14.05 


210.6^. 

213.7^ 


13.06* 

12.59^ 



a at B c . m . = 230 MeV. 
b at E c . m . = 300 MeV. 



light and medium mass systems as well as highly mass- 
asymmetric systems fusion generally occurs immediately 
above the ion-ion potential barrier, while in heavier sys- 
tems there is an energy range above the barrier where 
capture does not occur. The energy difference between 



the DC-TDHF potential barrier and the interaction bar- 
rier is the extrapush energy introduced by Swiatecki in 
a macroscopic model [36[; in addition to the work pre- 
sented here, this phenomenon has recently been studied 
for heavy and nearly symmetric reaction partners using 
the TDHF method E3- 



B. Dynamic excitation energy E*(R) 

In this subsection, we examine the dynamic excitation 
energy E*(R(t)), computed according to Eq. (gj), during 
the initial stages of the collision. Of particular interest is 
the excitation energy at the capture point, E*, which will 
influence the outcome of the reaction. In Fig. [5] the pre- 
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FIG. 5. (color online) Pre-compound excitation energy as 
a function of the internuclear distance R, calculated at zero 
impact parameter with the DC-TDHF method at four _E c . m . 
energies. Compared are the neutron-rich system 132 Sn+ 96 Zr 
(solid lines) and the stable system 124 Sn+ 96 Zr (dashed lines). 



compound excitation energies are shown as a function of 
the internuclear distance R; this represents our first mi- 
croscopic calculation for neutron-rich systems. When the 
two nuclei are far apart, the excitation energy is zero (this 
provides a good test for the numerical accuracy of the 
DC-TDHF calculation). As the two ions approach each 
other the excitation energy increases rapidly and reaches 
values between 30 - 90 MeV for the given range of cm. 
energies. It is interesting to note that at -E c .m. = 220 
MeV TDHF theory predicts that the two ions bounce 
back despite the fact that they are almost 10 MeV above 
the corresponding potential barrier, i.e. at this energy we 
have predominantly deep-inelastic reactions rather than 
capture. This is due to the fact that a large part of the 
incoming cm. energy was converted to internal excita- 
tion E* such that the collective energy does not suffice 
any more to surmount the barrier. This feature was also 



shown in the figure for the corresponding heavy-ion in- 
teraction potential (dotted line in Fig. [3]). 

In Fig.|5]we show the excitation energy E* at the cap- 
ture point as a function of E c _ m , . The capture point is de- 
fined as the distance R inside the barrier region where the 
collective kinetic energy, Eq. ([5]), becomes zero. Because 
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FIG. 6. (color online) Pre-compound excitation energy at the 
capture point, E*, as function of the cm. energy, as predicted 
by DC-TDHF for a central collision. 

of the very elongated shape of the dinuclear composite 
system at the capture point, E* is systematically lower 
than one would expect for the compound nuclei 228 > 220 Th 
in their ground state (dashed line in Fig. [6]) 

E* = E c . m . + Q gg . (10) 

The last expression can be derived from reaction kine- 
matics. The Q-values are obtained from measured bind- 
ing energies of the reaction partners and the compound 
nuclei. For the two systems considered here the Q-values 
are almost identical: Q gg = -188.7 MeV for 132 Sn+ 96 Zr 
and Q gg = -188.4 MeV for 124 Sn+ 96 Zr, hence we have 
drawn only one curve for both systems (dashed line in 
Fig.[B]). We observe that the excitation energy E* at the 
capture point is somewhat lower for the 124 Sn+ 96 Zr sys- 
tem in the energy range E c . m . = 230 — 250 MeV; at higher 
£-c.m. energies, their excitation energies are almost identi- 
cal. In this context, we would like to mention that recent 
microscopic calculations [38| have shown that the tem- 
perature (excitation energy) of the actinide compound 
nuclei will strongly influence the height of their fission 
barriers. 



Capture and deep-inelastic cross section for 

132 Sn+ 96 Zr 



Previously, we have studied heavy-ion fusion of the 
neutron-rich system 132 Sn+ 64 Ni using the DC-TDHF 



method [2l|, |22|. In that case, the fission barrier of the 
compound system is so high that the fission contribution 
is negligible at energies near the Coulomb barrier. By 
contrast, the compound nuclei for the systems studied 
in the present paper, 132,124 Sn+ 96 Zr, are the actinides 
228,220 Th with a fission barrier of only about 6 MeV. 
We therefore expect sizable fission competition, and the 
evaporation residue cross section will be quite small. De- 
pending upon beam energy and impact parameter, the 
dominant reaction channels are deep-inelastic and cap- 
ture reactions. In general, central collisions and collisions 
with relatively small impact parameter result in capture 
(one fragment in the exit channel), while at larger im- 
pact parameters the system tends to disintegrate into 
two fragments after some mass and charge transfer (deep- 
inelastic reactions). Regarding the capture channel, the 
composite system will eventually decay by nucleon and 
photon emission or by fission. This long-time evolution of 
the composite system is beyond the scope of TDHF due 
to the absence of quantum decay processes and transi- 
tions. In Fig. [7] we show total capture and deep-inelastic 
cross sections for the neutron-rich system 132 Sn+ 96 Zr. 
Let us first discuss the results obtained with the DC- 



gies -Ec.m. < 230 MeV correspond to a dinuclear system 
where both nuclei maintain separate cores, i.e. deep- 
inelastic reactions. The DC-TDHF capture cross section 
at i?c.m. — 230 MeV has been calculated from the energy- 
dependent heavy-ion potential at 230 MeV. Similar calcu- 
lations were carried out at energies -E c .m. = 240 and 250 
MeV. These capture cross sections are marked by filled 
triangles, and the solid line represents an interpolation 
between the data points. No restrictions were applied to 
the sum of partial waves L in Eq. ([3]), and this curve is 
therefore marked "all L" . The capture cross sections at 
energies below 230 MeV were obtained from the heavy- 
ion potential at 230 MeV because this is the lowest po- 
tential barrier that leads to capture. The deep-inelastic 
cross section at E CAn , = 220 MeV was calculated using 
the energy-dependent heavy-ion potential at 220 MeV. A 
similar calculation was carried out at energy -E c . m . = 210 
MeV. These deep-inelastic cross sections are marked by 
open triangles, and the dashed line represents an inter- 
polation between these data points. The capture cross 
sections at energies below 210 MeV were calculated from 
the heavy-ion potential at 210 MeV because this is the 
lowest potential barrier predicted by DC-TDHF. 
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FIG. 7. (color online) Total cross section for capture (solid 
line) and for deep-inelastic reactions (dashed line) for the 
neutron-rich system 132 Sn+ 96 Zr calculated with the DC- 
TDHF method as function of E c . m .. Total capture cross sec- 
tions predicted by unrestricted TDHF calculations are also 
given (red dots with error bars). For details, see the text. 



TDHF method, which can be used at energies -E c .m. be- 
low and above the potential barriers. From a compari- 
son of the heavy- ion potentials in Fig. 3] with the nuclear 
density distributions at the distance of closest approach, 
see Fig. [TJ we conclude that only potentials calculated 
at energies B c . m . > 230 MeV lead to a true composite 
system with overlapping cores, i.e. a capture reaction. 
By contrast, the heavy-ion potentials calculated at ener- 



FIG. 8. Total cross section for capture (solid line) and for 
deep-inelastic reactions (dashed line) for the neutron-rich sys- 
tem 132 Sn+ 96 Zr calculated with the DC-TDHF method as 
function of -E c .m. ■ The experimental capture-fission cross sec- 
tions are taken from Ref. [j|. 

At energies above the potential barriers (no barrier 
tunneling), we have also carried out unrestricted TDHF 
runs with impact parameters in the range of — 4 fm. 
By examining the density contours as a function of time, 
one can easily distinguish between capture events (one 
fragment in the exit channel) and deep-inelastic reac- 
tions (two fragments). At E c , m , = 250 MeV we find 
that impact parameters 6 = — 3.25 fm result in cap- 
ture, while impact parameters b > 3.50 fm lead to deep- 
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inelastic reactions. Using the sharp cut-off model, the 
capture cross section is given by <r cap t = '^b\ lax w hh 
bmax = 3.375 ± 0.125 fm. This cross section is shown 
by a red dot in Fig. [71 with the corresponding theoretical 
error bar arising from the impact parameter spacing. The 
impact parameter b max = 3.375 fm corresponds to an or- 
bital angular momentum quantum number L max = 87. 
If we use this angular momentum cut-off in the IWBC 
method, we obtain the curve labeled L max = 87; as 
we can see, both methods yield the same capture cross 
section. We have carried out a similar calculation at 
E c , m . = 230 MeV which yields b max = 2.05 ±0.15 fm and 
L max = 51. The corresponding TDHF and DC-TDHF 
capture cross sections are also plotted in Fig. [7] 
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FIG. 9. (color online) Total cross section for capture (solid 
line) and for deep- inelastic reactions (dashed line) for the sta- 
ble system 124 Sn+ 96 Zr calculated with the DC-TDHF method 
as function of E c , m .. The red dot represents the total capture 
cross section at E c , m . — 250 MeV predicted by an unrestricted 
TDHF run. For details, see text. 



In Fig. [5] we compare the DC-TDHF cross sections 
for deep-inelastic and capture reactions to experimen- 
tal capture-fission cross sections measured at HRIBF 
with a radioactive 132 Sn beam Q. Because the fission 
probability Pa ss < 1, our calculated capture cross sec- 
tions should be regarded as an upper limit for the mea- 
sured capture-fission cross sections. According to our 
unrestricted TDHF calculations, the dominant reaction 
channels at energies -E c .m. < 230 MeV are the deep- 
inelastic and quasi-elastic channels. In fact, our den- 
sity plots in Fig. [TJ reveal that at the lowest measured 
energy E c , m , — 195 MeV the nuclear surfaces barely 
touch. Any fission from such an event would have to 
arise from sub-barrier neutron-transfer and should be 
negligible compared to capture fission at higher ener- 
gies. We therefore make the conjecture that the bulk 
of the low-energy experimental data in Fig. [5] represent 



deep-inelastic and quasi-elastic events masquerading as 
capture-fission. Indeed, because of the limited mass res- 
olution in the HRIBF experiments Q it has been diffi- 
cult to separate the DIC component from capture-fission. 
Further experiments with an improved fission fragment 
detector are planned |39j . 



D. Capture and deep-inelastic cross section for 

124 Sn+ 96 Zr 



In Fig. [9] we examine the properties of the stable system 
124 Sn+ 96 Zr. Plotted are total cross sections for capture 
and for deep-inelastic reactions calculated with the DC- 
TDHF method as function of the cm. energy. Like in the 
corresponding neutron-rich system, only the heavy-ion 
potentials calculated at energies E CAn , > 230 MeV lead to 
capture while potentials at lower energies are associated 
with deep-inelastic channels. We have also carried out 
unrestricted TDHF runs for this system at -E c .m. = 250 
MeV. At impact parameter b = 3.5 fm we find a deep- 
inelastic reaction, and at b = 3.0 fm we obtain capture. 
Again, using the sharp cut-off model with b max = 3.0 fm 
we obtain the cross section shown by a red dot in Fig. ® 
This impact parameter corresponds to an orbital angular 
momentum quantum number L max — 77. If we use this 
angular momentum cut-off in the DC-TDHF method, we 
obtain the curve labeled L max = 77; as we can see, both 
methods yield the same capture cross section. 

1000 r — i 1 1 1 1 1 r 



124 Sn + 96 Zr 



H 100: 



10 L 



i 




f L =77 

max 

f 

f 



DC-TDHF: capture 

- - DC-TDHF: deep-inel. 
• exp: capture -fission 



200 



220 



240 



260 



E (MeV) 

cm. v 7 



FIG. 10. Total cross section for capture (solid line) and for 
deep-inelastic reactions (dashed line) for the stable system 
124 Sn+ 96 Zr calculated with the DC-TDHF method as func- 
tion of Ec.m.. The experimental capture-fission cross sections 
are taken from Ref. [jj]. 

In Fig. [TU] we compare the DC-TDHF cross sections 
for deep-inelastic and capture reactions to experimental 
capture-fission cross sections measured at HRIBF. Again, 



our calculated capture cross sections represent an upper 
limit for the measured capture-fission data. The agree- 
ment between theory and experiment is quite remarkable 
in view of the fact that we employ a fully microscopic 
theory based on a given energy functional, with no ad- 
justable parameters related to the capture process. 



IV. CONCLUSIONS 

In this paper we have studied deep-inelastic and cap- 
ture reactions for the neutron-rich system 132 Sn+ 96 Zr at 
energies in the vicinity of the barrier. This is by far 
the heaviest neutron-rich system we have investigated 
using both unrestricted TDHF and DC-TDHF methods. 
To elucidate any special properties neutron-rich systems 
might possess, we have compared a number of observables 
to those of the stable system 124 Sn+ 96 Zr. The dynamic 
microscopic calculations are carried out on a 3D Carte- 
sian lattice, and they require a large amount of CPU 
time, particularly with the added DC method. 

A contour plot of the mass density of the dinuclear sys- 
tem shows clearly a transition from two separate cores at 
lower energies to a shape configuration with overlapping 
cores or a single-core at energies E CAn , > 230 MeV. A 
study of the dynamic quadrupole moment Q2o{t) shows 
that even at E c , m , = 300 MeV, the intrinsic quadrupole 
moment is 3 times larger than that of the deformed com- 
pound nucleus 228 Th during the initial stages of the col- 
lision. We also calculate the heavy-ion interaction poten- 
tial V(R), and we demonstrate that in these very heavy 
systems the barrier height and width increase dramati- 
cally with increasing beam energy E cra . We find that the 
potential barriers of the neutron-rich system 132 Sn+ 96 Zr 
are systematically 1 — 2 MeV higher than those of the 



stable system. By contrast, we observe large differences 
(9 MeV) in the interaction barriers of the two systems 
which can be deduced from unrestricted TDHF runs. We 
then examine the dynamic excitation energy E*(t) dur- 
ing the initial stages of the collision and compare it to 
the expression E* = E CnTn _ + Q gg (deduced from reaction 
kinematics) which assumes that the combined system is 
in its ground state. Finally, capture cross sections for the 
two reactions are analyzed in terms of dynamic effects, 
and a comparison with recently measured capture-fission 
data is given. 

One of the major open questions in the reactions of 
neutron-rich nuclei is the dependence of the barrier on 
isospin T z = (Z — N)/2. To reveal possible systematic 
trends requires additional theoretical and experimental 
studies with a wide variety of projectile and target com- 
binations which are expected to become available at cur- 
rent and future RIB facilities. To be able to pin down 
the isospin dependence in a fully microscopic theory, it 
is probably best to choose collision partners which are as 
simple as possible: projectile and target nuclei should be 
spherical in their ground state, and the compound nu- 
cleus should have a high fission barrier so that the fission 
component can be ignored (at least at lower beam ener- 
gies) . A desirable reaction system of this kind appears to 
be 132 Sn+ 40 ' 48 ' 54 Ca, and we are planning to investigate 
these reactions in the future. 



ACKNOWLEDGMENTS 

This work has been supported by the U.S. Department 
of Energy under Grant No. DE-FG02-96ER40963 with 
Vanderbilt University, and by the German BMBF under 
Contracts No. 06F131 and No. 06ER142D. 



[1] J. F. Liang et al, Phys. Rev. Lett. 91, 152701 (2003); [12] 

96, 029903(E) (2006). 
[2] J. F. Liang. D. Shapira, C. J. Gross, R. L. Varner, J. R. 

Beene, P. E. Mueller, and D. W. Stracener, Phys. Rev. 

C 78, 047601 (2008). [13] 

[3] A. M. Vinodkumar et al, Phys. Rev. C 78, 054608 

(2008). 
[4] S. Hofmann et al, Eur. Phys. J. A 14, 147 (2002). [14] 

[5] Yu. Ts. Oganessian et al, Phys. Rev. Lett. 104, 142502 

(2010). [15] 

[6] T. Ichikawa, A. Iwamoto, P. Moller, and A.J. Sierk, Phys. 

Rev. C 71, 044608 (2005). [16] 

[7] V. Yu. Denisov and W. Norenberg, Eur. Phys. J. A 15, 

375 (2002). [17] 

[8] G. Fazio, et. al, Eur. Phys. J. A 19, 89 (2004). 
[9] A. K. Nasirov, G. Giardina, G Mandaglio, M. Manga- [18] 

naro, F. Hanappe, S. Heinz, S. Hofmann, A. I. Muminov, 

W. Scheid, Phys. Rev. C 79, 024606-10 (2009). [19] 

[10] G. G. Adamian, N. V. Antonenko , and W. Scheid, Eur. 

Phys. J. A 41, 235 (2009). [20] 

[11] J. W. Negele, Rev. Mod. Phys. 54, 913 (1982). 



K. T. R. Davies, K. R. Sandhya Devi, S. E. Koonin, and 

M. R. Strayer, page 3 in "Treatise on Heavy-Ion Physics, 

Vol. 3 Compound System Phenomena" edited by D. A. 

Bromley, Plenum Press, New York, 1985. 

E. Chabanat, P. Bonche, P. Haensel, J. Meyer and R. 

Schaeffer, Nucl. Phys. A635, 231 (1998); A643, 441(E) 

(1998). 

M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. 

Phys. 75, 121 (2003). 

A. S. Umar and V. E. Oberacker, Phys. Rev. C 73, 

054607 (2006). 

Lu Guo, J. A. Maruhn, P.-G. Reinhard, and Y. 

Hashimoto, Phys. Rev. C 77, 041301(R) (2008). 

Kouhei Washiyama and Denis Lacroix, Phys. Rev. C 78, 

024610 (2008). 

David J. Kedziora and Cedric Simenel, Phys. Rev. C 81, 

044613 (2010). 

Y. M. Engel, D. M. Brink, K. Goeke, S. J. Krieger, and 

D. Vautherin, Nucl. Phys. A 249, 215 (1975). 

A. S. Umar and V. E. Oberacker, Phys. Rev. C 74, 

02160RR) (2006). 



10 



[21] A. S. Umar and V. E. Oberacker, Phys. Rev. C 74, 

061601(R) (2006). 
[22] A. S. Umar and V. E. Oberacker, Phys. Rev. C 76, 

014614 (2007). 
[23] A. S. Umar and V. E. Oberacker, Phys. Rev. C 77, 

064605 (2008). 
[24] A. S. Umar and V. E. Oberacker, Eur. Phys. J. A 39, 

243 (2009). 
[25] A. S. Umar, J. A. Maruhn, N. Itagaki, and V.E. Ober- 
acker, Phys. Rev. Lett. 104, 212503 (2010). 
[26] A. S. Umar, V. E. Oberacker, J. A. Maruhn, and P.-G. 

Reinhard, Phys. Rev. C 81 064607 (2010). 
[27] A. S. Umar and V. E. Oberacker, Phys. Rev. C 74, 

061601(R) (2006). 
[28] G. H. Rawitscher, Phys. Rev. 135, 605 (1964). 
[29] K. Hagino and Y. Watanabe, Phys. Rev. C 76, 021601 (R) 

(2007). 



[30] R. Y. Cusson, P.-G. Reinhard, M. R. Strayer, J. A. 

Maruhn, and W. Greiner, Z. Phys. A 320, 475 (1985). 
[31] A. S. Umar, V. E. Oberacker, J. A. Maruhn, and P.-G. 

Reinhard, Phys. Rev. C 80, 041601(R) (2009). 
[32] A. S. Umar and V. E. Oberacker, J. Phys. G: Nucl. Part. 

Phys. 36, 025101 (2009). 
[33] G. R. Satchler and W. G. Love, Phys. Rep. 55, 183 

(1979). 
[34] M. J. Rhoades-Brown, V. E. Oberacker, M. Seiwert, and 

W. Greiner, Z. Phys. A 310, 287 (1983). 
[35] G. Bertsch, J. Borysowicz, H. McManus, and W. G. Love, 

Nucl. Phys. A 284, 399 (1977). 
[36] W. Swiatecki, Nucl. Phys. A 376, 275 (1982). 
[37] Cedric Simenel, Benoit Avez, Cedric Golabek, 

( [http://arxiv.org/abs/0904.2653vl | ) . 
[38] J. C. Pei, W. Nazarewicz, J. A. Sheikh, and A. K. Ker- 

man, Phys. Rev. Lett. 102, 192501 (2009). 
[39] D. Shapira, ORNL, private communication. 



